AD-fll44  876  A  PERTURBATION  SOLUTION  FOR  RIGID  DUMBEEll  SUSPENSIONS.  1/1 
IN  STEADV  SHEAR  FLOW(U)  WISCONSIN  UNIV-HADISON 
MATHEMATICS  RESEARCH  CENTER  S  KIM  ET  AL.  JUN  84 
UNCLASSIFIED  MRC-TSR-2700  DAAG29-80-C-0041  F/G  12/1  NL 


line  FILE  COEi 


«.**.■  1  .■  \W9.9  v  • w  *115 1  «r  v*.  j'j f.i  v 


■  .»  v  V  \> V1^ 


Mathematics  Research  Center 
*  University  of  Wisconsin— Madison 
610  Walnut  Street 
Madison,  Wisconsin  53705 


June  1984 


(Received  January  4,  1984) 


Approved  for  public  release 
Distribution  unlimited 


OTIC. 

electe 

AU8  2  9 1984 


m 


Sponsored  by 

U.  S.  Army  Research  Office 
P.  0.  Box  12211 
Research  Triangle  Park 
North  Carolina  27709 


National  Science  Foundation 
Washington,  D.  C.  20550 


08  24  047, 

A-V-.'- A  A  A  -VV-V- A.'-  .*•  AA  -N  A  A 


_  • 


.'•A 

A'A'.V1 

>7$ 


UNIVERSITY  OF  WISCONSIN-MADISON 
MATHEMATICS  RESEARCH  CENTER 

A  PERTURBATION  SOLUTION  FOR  RIGID  DUMBBELL 
SUSPENSIONS  IN  STEADY  SHEAR  FLOW 

Sangtae  Kim*'^  and  Xijun  Fan**'^ 

Technical  Summary  Report  #2700 
June  1984 

ABSTRACT 

The  orientation  distribution  function  of  a  dumbbell  undergoing  rotary 
Brownian  motion  (with  time  constant  X)  in  steady  shear  flow,  and  the 
viscosity  and  normal  stress  coefficients  are  computed  by  a  perturbation 
expansion  in  small,  dimensionless  shear  rates  Xy.  It  is  found  that  the 
series  expansions  for  both  the  viscosity  and  the  first  normal  stress 
coefficient  have  identical  radii  of  convergence  at  Xy  =  0.81.  By  using 
analytical  continuation,  the  solution  was  extended  to  Xy  ■  1.5,  and  checked 
with  established  results.  Our  results  agree  with  earlier  perturbation 
calculations  of  Kirkwood  and  Plock  (1956),  Bird  and  Warner  (1971)  and 
numerical  calculations  of  Stewart  and  Sorenson  (1972). 
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SIGNIFICANCE  AND  EXPLANATION 


Polymeric  solutions  play  a  central  role  in  the  manufacture  of  many 
plastics  and  especially,  artificial  fibres.  A  large  research  effort  over 
decades  has  therefore  been  directed  towards  the  prediction  of  the  mechanical 
properties  of  such  solutions  which  are  important  in  manufacture,  in 
particular,  their  viscosity  and  their  normal  stress. 

It  has  turned  out  that  many  features  found  in  a  plot  of  the  viscosity  and 
normal  stress  coefficients  vs.  shear  rate  for  polymeric  solutions  can  be 
modeled  qualitatively  by  a  suspension  of  rigid  dumbbells.  Thus  this  model  has 
been  quite  popular.  Despite  many  attempts,  exact  analytical  solutions  have 
never  been  found  for  this  problem  although  numerical  and  asymptotic  results 
are  available. 

The  following  investigation  establishes  that  the  asymptotic  solution  for 

weak  flows  (or  strong  Brownian  diffusion)  is  accurate  only  for  dimensionless 

•  • 

shear  rates  Xy  less  than  0.81  where  Y  is  the  shear  rate  and  X  is  a 
material  time  constant  for  the  rigid  dumbbell  model. 
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INTRODUCTION 


The  rigid  dumbbell  in  a  Newtonian  solvent  under  hydrodynamic  and 
?-cwnian  forces  is  a  simple  model  of  solutions  of  rigid  macromolecules  with 
both  theoretical  and  pedagogical  applications.1 In  this  note,  we  show  that 

S' 

expansion  solutions  in  small,  dimensionless  shear  rates, (Ay,  have  a  finite 

✓ 

radius  of  convergence.  We  also  show  that  these  solutions  can  be  analytically 
continued  to  a  much  larger  domain  by  standard  techniques  for  accelerating  the 
convergence  of  series. 

V  . 

STATEMENT  OF  THE  PROBLEM 

Consider  a  dilute  suspension  (number  of  particles  per  unit  volume, 
n  <<  1)  of  rigid  axi symmetric  particles,  where  hydrodynamic  interaction  is 
negligible  between  particles.  The  material  functions  can  be  determined  from 
the  orientation  distribution  function,  ip(6,f),  where  *(6,$)  ded$  gives  the 
fraction  of  particles  with  the  axis  oriented  between  (e,  e  +  de)  and 
(4>,  i>  +  d4>).  0  and  4>  are  the  usual  spherical  azimuthal  and  longitudinal 

angles.  The  equation  of  continuity  and  the  equations  of  motion  combine  to 
give  the  "diffusion  equation"  for  ip,  which  for  steady  shear  flow  becomes^ 

A(sTnF)  "  6A^fi(sTnF)  *  0  *  (1) 


with  the  operators  A,  and  fi  defined  as 
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We  have  defined  our  coordinate  system  such  that - =  -y  is  the  only  non- 

ay 

vanishing  element  of  the  velocity  gradient  tensor  Vv_.  For  rigid  dumbbells  of 
length  L,  and  bead  friction  coefficient  the  rotational  time  constant  is 
given  by  X  »  cL2/(12kT).  The  expression  for  X  for  other  axi symmetric  shapes 
is  given  by  Brenner.2 

A  general  analytic  solution  of  (1)  has  never  been  found  although  exact 

A 

solutions  of  related  equations  such  as  Boeder's  equation  are  known. 

Our  objective  is  to  solve  equation  (11  for  ^  using  a  perturbation  series 
in  Xy  and  to  use  the  results  to  calculate  the  viscosity  n  and  the  first 
normal -stress  coefficient  ^  vanishes  for  this  rigid  dumbbell  model): 

n  -  ns  *  I  nkTx<2(P°  -  P°)  -  ?\  cos  2*>+  (2) 


T.  «  Hiili  <P2  sin  24>> 
y 


THE  PERTURBATION  SOLUTION 

We  expand  ip  in  powers  of  \y 


BO 

*(©,-♦)  =  i-  I  (6X^)k  4>k(  6 , 4>) 
*  k*0  K 


and  substitute  into  the  diffusion  equation  (1)  to  get  the  following  hierarchy 
of  equations  for  the 
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(8) 


- -■  I  [A(l ,1 ,2k)  -  I  A(l,2,2k)  -  £  A(2,2,2k)  ](6x-*)2k 

nllTX  k*0  6  5 


(-l)k  ak(X^)2k 
k*0  * 


— - - «  I  12A(2,2,2k+l)(6Xt)2k 

6nkTX2  k*° 

(9) 

5  I  (-l)k  bk(Xt)2k 
k-0  K 

The  first  twenty  coefficients,  a*  and  b^,  are  given  in  Table  1.  The  values 
for  aQ,  aj^  and  a2  agree  with  Kirkwood  and  Plock5  and  Bird  and  Warner.6 

DISCUSSION 

We  note  that  the  ratios  aic+j/a^  and  b^+i/b^  approach  an  asymptotic  value 
with  large  k,  i.e.  approximately  1.521006.  This  result  was  tested  by 
generating  a^.and  b|c  up  to  k  *  85,  and  it  was  found  to  hold  for  these  higher 
coefficients  as  well.  Therefore,  the  coefficients  approach  asymptotically  to 
the  coefficients  of  a  geometric  series  in  (Xy)2  with  a  radius  of  convergence 
of  ( 1 . 521006 ) ~ 1 /2  s  o.81.  Therefore,  our  series  solution  converges  only  for 
Xy  <  0.81.  However,  since  the  coefficients  approach  those  of  a  geometric 
series,  we  tried  analytical  continuation. 

First,  we  replace  the  higher  coefficients  a^  and  bfc  for  k  >  N+l  with 
their  asymptotic  counterparts  a^  and  b|<: 


Table  1 


Coefficients  for  the  Series  Solutions 


0.5142  8571 
0.6888  3117 
1.0173  986 
1.5354  447 
2.3303  602 

0.3542  3678*101 
0.5387  1583*10! 
0.8193  7480*101 
0.1246  2985*102 
0.1895  6859*102 

0.2883  4481*102 
0.4385  8958*102 
0.6671  2099*102 
0.1014  7310*103 
0.1543  4667 *103 

0.2347  7055*103 
0.3571  0009 *103 
0.5431  7067  *103 
0.8261  9517  *103 
0.1256  6924*104 


1.0857  143 
1.5475  724 
2.3156  606 
3.5063  866 
5.3265  547 

0.8098  9823*101 
0.1231  7694*102 
0.1873  5352*102 
0.2849  7327*102 
0.4334  6025*102 

0.6593  1852*102 
0.1002  8628*103 
0.1525  4144*103 
0.2320  2467*103 
0.3529  2344*103 

0.5368  1773*103 
0.8165  3198*103 
0.1241  9941 *104 
0.1889  1476*104 
0.2873  5068*104 
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4>i  ■  —  P;  sin  24 

i  12  c 
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k  •  2,  3,  4  ... 


As  pointed  out  by  Kirkwood  and  Plock®  and  by  Bird  and  Warner®  the  4^  may  be 
expressed  as  a  linear  combination  of  the  spherical  harmonics,  due  to 
the  fact  that 


MpMH  It)  •  -n(n+l)  ;♦ 

n  ‘cos  m®  n 'cos  md> 


m*.  "t2  %2  nj.j.cos  l 

j.I-2  k-L  WkL,1„ 


mj 

The  an(t  are  tabulated  in  Table  1  of  reference  (6)  and  Table  11.4-1  of 
reference  (2).  Furthermore,  from  (7)  it  follows  that  cos  m$>  terms  appear 
only  in  the  ^  with  k  even,  while  sin  m$  terms  appear  only  with  k  odd. 
Therefore,  the  material  functions  given  by  equations  (2)  and  (3)  reduce 
to  a  power  series  in  (Xy)2  because  of  the  orthogonality  with  respect  to  the 
4  integration.  The  perturbation  results  up  to  ^  are  given  in  reference  (2) 
We  have  programmed  this  procedure  (i.e.  equations  4,  6  and  7),  by 
generating  a  recursion  formula  for  A(m,n,k),  where  the  m,n,k  element  is  the 
coefficient  of  ^(n-l)  ^or  *k’  "^e  efficients  generated  by  the  program 

for  4>j,  $2  and  $3  matched  those  given  in  reference  (2).  Because  of 
orthogonality,  only  a  few  terms  make  a  contribution  to  the  material 
functions,  i.e.: 


{-l)k  ak(xt)2k  +  2  (-^k  V>t)J 


k«N+l 


S  (-l)k  bk(XT)2k  +  2  M>k  Mx*): 


6nkTX 


k*N+l 


with  &k  ■  c1p”2k"2,  Cj  *  0.1879, 


p  *  0.8108., 


bk  «  C2p“Zk”2,  c2  «  0.3582... 


p  *  0.8108... 


In  the  key  step,  the  series  Is  replaced  by  the  equivalent  continuation 


functions  so  that 


2  (-l)k[ak-ak](xt)2k  +  Cj[p2  +  (x-r)2]' 


2  (-l)k[bk-bk](xt)Zk  +  c2[pz  +  (xt)z]‘ 


6nkTX  k*0 


By  comparing  with  the  exact  numerical  (collocation)  solutions  of  Stewart  and 
Sorensen7  we  found  that  N  =  10  was  the  optimum  value.  However,  N  *  8,  9,  11, 


and  12  also  gave  results  which  were  accurate  to  four  significant  figures  for 
Xy  <  1.  Table  2  shows  the  comparison  between  their  collocation  solution  and 


~  — . . <  ^  ^ 


Comparison  of  Stress  Functions 


0 

0.1 

0.125 

0.25 

0.3333 

0.5 

0.75 

1.0 

1.5 

2.0833 


(n-ns)/nkTX  ^/nkT^ 

equation  (12)  exact  equation  (13)  exact 

—  A  a  10  solution8  hi  «  10  solution8 


1 

1 

6/5 

6/5 

0.9949 

0.9949 

1.1872 

1.1872 

0.9921 

0.9921 

1.1801 

1.1801 

0.9703 

0.9703 

1.1252 

1.1252 

0.9502 

0.9502 

1.0749 

1.0749 

0.9029 

0.9029 

0.9588 

0.9588 

0.8306 

0.8306 

0.7874 

0.7874 

0.7676 

0.7676 

0.6467 

0.6467 

0.6604 

0.6735 

0.4596 

0.4565 

9.2 

0.5998 

-5.5 

0.3289 

8Col location  calculations  of  Stewart  and  Sorensen.7 


our  extended  solution  (with  N  *  10,  c^  »  0.18799972,  *  0.35822745, 

\ 

and  p  *  0.81082443). 

ACKNOWLEDGEMENTS:  The  authors  thank  Professor  R.  B.  Bird  for  bringing  our 
attention  to  this  oroblem  and  Professor  W.  E.  Stewart  and  Dr.  J.  P.  Sorensen 
for  furnishing  details  of  their  numerical  calculations. 


-8- 


RE- “R^CES 


(1)  w.  Kuhn,  Kolloid.  Z.,  68,  2-15  (1934). 

(2)  R.  B.  Bird,  0.  Hassager,  R.  C.  Armstrong,  and  C.  F.  Curtiss. 

Dynamics  of  Polymeric  Liquids:  Volume  2.  Wiley,  New  York  (1977). 

(3)  H.  Brenner,  Int.  J.  of  Multiphase  Flow,  1_,  195  (1974). 

(4)  H.  Brenner,  J.  Coll.  Int.  Sci.,  34,  103  (1970). 

(5)  J.  G.  Kirkwood  and  R.  J.  Flock,  J.  Chem.  Phys.,  24,  665  (1956). 

(6)  A.  B.  Bird  and  H.  R.  Warner,  Jr.,  Trans.  Soc.  Rheol . ,  15,  741  (1971) 

(7)  W.  E.  Stewart  and  J.  P.  Sorensen,  Trans.  Soc.  Rheol.,  16,  1  (1972). 


t 


SECURITY  CLASSIFICATION  OP  THIS  PAGE  (Whm I  dm*  Knltrod) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


.  RECIPIENT  $ CATALOO  NUMBER 


4.  TITLE  (ond  Submit) 

A  PERTURBATION  SOLUTION  FOR  RIGID  DUMBBELL 
SUSPENSIONS  IN  STEADY  SHEAR  FLOW 


.  author/*; 

Sangtae  Kin  and  Xijun  Fan 


PERPORMINO  ORGANIZATION  NAME  AND  AODRCM 

Mathematics  Research  Center,  University  of 
610  Walnut  Street  Wisconsin 

Madison.  Wisconsin  53706 _ 


II.  CONTNOLLINO  OFFICE  NAME  AND  AODRCM 

See  Item  18  below. 


ADDRESS/!/  4 ft* ran  I  I 


Centrolltnd  OlHeo) 


4.  PERPORMINO  ORO.  REPORT  NUMBER 


•.  CONTRACT  OR  orant  number/*; 

DAAG29-80-C-00  41 
CPE-8104705 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  A  WORK  UNIT  NUMBERS 

Work  Unit  Number  2  - 
Physical  Mathematics 


I*.  REPORT  DATE 

June  1984 


IS.  number  op  paoes 

9 


IS.  SECURITY  CLASS,  (ol  tklo  report; 


UNCLASSIFIED 


DECL  ASSI  PI  C  ATI  ON/  DOWNOR  ADINO 
SCHEDULE 


S.  DISTRIBUTION  STATEMENT  /•/  mi*  Report) 

Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (ol  »•  obotroct  onlocod  In  Block  to,  II  dftforont  Irom  Report; 


IS.  SUPPLEMENTARY  NOTES 

U.  S.  Army  Research  Office  National  Science  Foundation 

P.  0.  Box  12211  Washington,  D.  C.  20550 

Research  Triangle  Park 

North  Carolina  27709  _ 


IS.  KEY  WORDS  (Conllnu*  on  rorotoo  old*  II  noooooot y  oN  tdontlfr  b y  block  mmbot) 

rigid  dumbbells 

rotary  Brownian  diffusion 

suspension 


20.  ABSTRACT  (Conllnu*  on  rororoo  old*  II  noeooooty  ond  Identity  by  block  mmbot) 

The  orientation  distribution  function  of  a  dumbbell  undergoing  rotary 
Brownian  motion  (with  time  constant  X)  in  steady  shear  flow,  and  the 
viscosity  and  normal  stress  coefficients  are  computed  by  a  perturbation 
expansion  in  small,  dimensionless  shear  rates  \y .  It  is  found  that  the 
series  expansions  for  both  the  viscosity  and  the  first  normal  stress 


PORM  in* 
I  JAN  72  H/J 


COITION  OP  I  NOV  SS  IS  OBSOLETE 


3 


UNCLASSIFIED _ 

SECURITY  CL  ASSI  PI  CAT  ION  OP  THIS  PAGE  (When  Dot*  Bntorod) 


».  *«  V  *.  _  m  ,1  »  *_l»  *  Jk  . 


20.  ABSTRACT  -  cont'd. 


coefficient  have  identical  radii  of  convergence  at  Xy  =  0.81.  By  using 
analytical  continuation,  the  solution  was  extended  to  Xy  =  1.5,  and  checked 
with  established  results.  Our  results  agree  with  earlier  perturbation 
calculations  of  Kirkwood  and  Plock  (1956)  ,  Bird  and  Warner  (1971)  and  numerical 
calculations  of  Stewart  and  Sorenson  (1972). 


